
NASA— TM-1 11253 


AIAA-95-1 292 

Euler/Navier-Stokes Flow Computations 
on Flexible Configurations for Stability 
Analysis 


G. Guruswamy and E. Tu 
NASA Ames Research Center 
Moffett Field, CA 


36th AIAA/ASME/ASCE/AHS/ASC 
Structures, Structural Dynamics 
and Materials Conference 
April 10-12, 1995/New Orleans, LA 

For permission to copy or republish, contact the American Institute of Aeronautics and Astronautics 
370 L'Enfant Promenade, S.W., Washington, D.C. 20024 




AIAA— 95— 1 292— CP 


EULER / N AVIER- STOKES FLOW COMPUTATIONS ON FLEXIBLE 
CONFIGURATIONS FOR STABILITY ANALYSIS 

Guru P. Guruswamy * and Eugene L. Tu** 

Applied Computational Aerodynamics Branch, MS 258-1 
NASA Ames Research Center, Moffett Field, California 94035-1000 


Abstract 

Longitudinal dynamic stability derivatives re- 
quired for design of aircraft are computed by us- 
ing the state-of-the-art numerical methods for wing- 
body configurations. The flow is modeled using the 
Euler/Navier-Stokes equations with turbulence models 
and solved using an efficient finite- difference scheme 
suitable for patched structured grids. Computations 
are made at a flow regime that is beyond the limits of 
the current linear methods mostly used for computing 
stability derivatives. Flow conditions include shock- 
waves and viscous dominated vortical flows. Effect of 
Mach number and angle-of-attack on stability deriva- 
tives are demonstrated for a typical wing-body config- 
uration. For the same configuration the effects of wing 
flexibility on the magnitude and phase angles of stabil- 
ity derivatives are also demonstrated. 

Introduction 

Stability derivatives are required for design of air- 
craft. Most of the current stability derivative com- 
putations are based on the linear aerodynamic theory 
which is not adequate for transonic/ vortical/separated 
flows where flows are dominated by strong nonlineari- 
ties, including moving shock waves. In addition to flow 
non-linearities, aircraft flexibility can significantly in- 
fluence stability derivatives. Earlier wind tunnel mea- 
surements made on the Boeing 2707 SST model show 
significant differences in stability derivatives between 
rigid and flexible configurations. 1 Those differences are 
more pronounced in the transonic regime. Typical re- 
sults in Ref. 1 show that the stability derivatives for 
the rigid and flexible configurations differ by as much 
as 30% at a transonic Mach number near 0.90. 

In order to accurately determine stability deriva- 
tives associated with strong flow non-linearities, current 
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practice is to conduct wind-tunnel tests. Wind-tunnel 
tests for transonic flows are not only expensive, but 
can also be time consuming. Furthermore, including 
important flexibility effects can make wind-tunnel ex- 
periments even more expensive and can often be more 
dangerous. As a result, parallel work using computa- 
tional methods is required to confirm and complement 
wind tunnel measurements. Computational models are 
considerably less expensive and can yield early useful 
results which can help wind tunnel tests. With this ef- 
fort, the overall cost of analysis and design of aircraft 
can be reduced considerably. 

In order to accurately model non-linear flows that 
involve moving shock waves, vortices and flow separa- 
tions, exact equations such as the Euler/Navier-Stokes 
equations are required. For the past two decades, Com- 
putational Fluid Dynamics (CFD) methods have made 
significant progress. Codes based on the Euler/Navier- 
Stokes equations coupled with structural equations are 
in routine use. At NASA Ames Research Center, 
significant effort has been made to develop a com- 
puter code, ENSAERO, that simultaneously solves the 
Euler/Navier-Stokes flow equations and modal/finite- 
element structural equations. 2 This code has been ex- 
tended to compute unsteady flow on complex rigid con- 
figurations such as wing-body-canard and wing-body- 
control configurations by using patched zonal grids. 3 ’ 4 
ENSAERO has a capability to compute aeroelastic- 
ity of wing-body configurations by using modal/finite- 
element structures with aeroelastically deforming built- 
in grids. 5 

In this work, the capability of ENSAERO to com- 
pute stability derivatives associated with complex flows 
is demonstrated for flexible configurations. Longitudi- 
nal dynamic stability derivatives required for design of 
aircraft are computed by using the state-of-the-art nu- 
merical methods for rigid and flexible wing-body con- 
figurations. Computations are made at a flow regime 
that is beyond the limits of the linear methods mostly 
used for computing the stability derivative. Flow con- 
ditions include shock- waves and viscous dominated vor- 
tical flows. Effect of Mach number and angle-of-attack 
on stability derivatives are demonstrated for a typical 
wing-body configuration. For the same configuration 
the effect of wing flexibility on the magnitude and phase 
angles of stability derivatives is also demonstrated. 
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Governing Aerodynamic Equations 

The strong conservation-law form of the Navier- 
Stokes equations are used for shock capturing purposes. 
The thin-layer version of the equations in generalized 
coordinates can be written as 6 

d r Q + d ( E + d v F + d c G = Re~ 1 d{S ( 1 ) 

where Q, E, F , G, and 5, are flux vectors in generalized 
coordinates. The following transformations are used in 
deriving Eq. (1). 

t — t 

£ = £(*,2 />M) 

T) = T](x,y,z,t) 

C = C(x,y,z,t ) 

It should be emphasized that the thin-layer approxi- 
mation is valid only for high Reynolds number flows 
and very large turbulent eddy viscosities invalidate the 
model. 

To solve Eq. (1), ENSAERO has time-accurate 
methods based on central-difference method. 7 

Moving Grids for Flexible Configurations 

One of the major difficulties in using the Navier- 
Stokes equations for computational aerodynamics lies 
in the area of grid generation. For steady flows, the 
advanced techniques such as zonal-grids are currently 
being used. However, grid-generation techniques for 
aeroelastic calculations, which involve moving compo- 
nents, are still in the early stages of development. 

In this work, an H-0 type grid topology is used. 
This grid topology is more suitable for general wing- 
body configuration. The base surface grid is gener- 
ated using the S3D code. 8 From the surface grid, the 
flowfield grid is generated using the elliptic grid code 
3DGRAPE. 9 For flexible configurations, the displace- 
ments at surface grid points are first computed. Then 
the flowfield grid is algebraically generated by redis- 
tributing the points along £ grid lines that are in the 
radial direction starting normal to the surface. The 
grid points are redistributed by moving them along 
C grid lines. Each grid point is displaced by a value 
that is equal to the surface displacement times the arc 
length measured from the outer boundary. Thus the 
flowfield grid will match with the structural displace- 
ments at the surface and will also remain stationary at 
the outer boundary. This concept was utilized in Ref, 
12 to quickly and efficiently generate the static grids at 
different canard deflections. 


Results 

Effect of Wing Flexibility on Time Responses 

In this paper a typical wing-body configuration 
shown in Fig. 1 is selected for demonstration. For this 
rigid configuration, the measured aerodynamic and sta- 
bility derivative data are given in Ref. 10 and 11. For 
this rigid configuration, steady, unsteady, and aeroelas- 
tic computations were made in Ref. 6 using ENSAERO. 

In Ref. 6, the wing-body configuration shown in 
Fig. 1 was modeled using a H-0 type grid and a grid 
refinement study was conducted. It was concluded that 
a grid size of 95 x 79 x 30 is adequate for computing 
pressures up to 12 deg angle of attack(AoA). The grid 
topology is shown in figure 1. 

Steady state computations which are required as 
initial conditions for unsteady computations were made 
at several flow conditions and results are compared with 
the experiment. Typical computations required about 
1000 steps for the residual to drop by three orders of 
magnitude. Figure 2 shows the comparison between 
the computed and the measured data at = 0.90, 
Re c = 1.52 x 10 6 , and a = 4.09 deg. The comparison 
is favorable. For this case, the flow is dominated by the 
presence of a moderate vortex and a weak shock-wave 
on the wing. This can be seen in the pressure-suction 
peaks near the leading edge in Fig. 2. 

In order to study the effect of flexibility on aerody- 
namic force responses, computations are made on rigid 
and flexible configurations. In this paper, effect of wing 
flexibility in time responses are studied. Computations 
are made for 0 to 12 deg AoA at Mach numbers varying 
from 0.90 to 1.10 for rigid and flexible configurations in 
which the wing is oscillating in its first torsional mode. 
From the previous studies in Ref. 12, it was found that 
at these conditions the flow is dominated by leading 
edge vortices. For both rigid and flexible computations, 
the complete configuration was oscillating in the pitch- 
ing mode with an amplitude of 1 deg about a pitch axes 
located at 2.45 root chords from the nose of the body. 
For the flexible case, the wing is oscillating in phase 
with the complete configuration with a tip amplitude 
of 1 deg. The reduced frequency, k = is °- 5 for 

both configuration and wing oscillations. 

At AoA = 0.0 computations were made using the 
Euler option of the code. The first grid spacing used 
is 0.005 of the root chord. The Euler computations 
required 2400 time steps per cycle. Two cycles were 
required during which the numerical transients disap- 
peared and a periodic solution was obtained. The lead- 
ing and trailing edge responses of the wing tip for the 
rigid and flexible configuration are shown in Fig. 3. 
The motion was started by first pitching downwards. 
As a result the tip is initially displaced upwards since 
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the axes of rotation for the configuration is forward of 
the wing-tip of leading edge. 

Figures 4 and 5 show the effects of flexibility on 
total lift and moment responses, respectively, for M = 
0.96. Both the lift and moment responses of the flexible 
case leads the rigid case in phase. The wing flexibility 
reduces magnitude of both the lift and moment. The 
flexible case has reduced the phase angle with respect 
to the motion when compared to the rigid case. Figures 
6 and 7 show the effect of Mach number on magnitude 
and phase angles of lift for both rigid and flexible cases, 
respectively. For both cases the magnitude increases 
and the phase angle decreases with the increase in Mach 
number. 

Computations similar to that carried out at AoA 
= 0 deg were also made at AoA = 4 deg by using the 
Navier-Stokes flow equations. A grid spacing of the or- 
der 0.00001 of root chord was used at the surface and 
the computations required 9000 time steps per cycle. 
Two cycles were required during which the numerical 
transients disappeared and a periodic solution was ob- 
tained. Figures 8 and 9 show the effect of flexibility on 
the magnitude and phase angles of the lift coefficient at 
various Mach Numbers ranging from 0.90 to 1.05. The 
magnitude of lift for both rigid and flexible cases de- 
creases with the increase in Mach number. The magni- 
tude of lift for the flexible case is always lower than that 
of the rigid case. Both phase angles in Fig. 9 increase 
with the increase in Mach number. This is opposite to 
what was observed at AoA = 0 deg. The flexible case 
again leads the rigid case for all Mach numbers, similar 
to the case of AoA = 0 deg. 

It was also of interest to study the effect of angle of 
attack on force responses. Computations were made at 
0, 4, 8 and 12 deg mean AoA. Figures 10 and 11 show 
the effect of AoA on magnitude and phase angles for 
both rigid and flexible cases, respectively. Lifts for both 
rigid and flexible cases increase until about 4 deg AoA 
and start decreasing. At all AoA, the rigid cases have 
higher magnitudes than the flexible cases. The differ- 
ence between rigid and flexible cases decreases slightly 
with the increase in AoA. Figure 11 shows that phase 
angles of the rigid cases increase until AoA of about 
4 deg and then start decreasing. For the flexible case 
the phase angle continuously decreases and the rate of 
decrease slightly increases after 4 deg AoA. Again the 
phase angle behaviors are influenced by vortical flows. 

The effect flexibility at various Mach numbers and 
AoA on the stability derivatives is studied in next sec- 
tion. 

Effect of Wing Flexibility on Dynamic Stability 

The dynamic stability parameters for the rigid and 
flexible wing-body configurations are predicted from 


the Fourier analysis presented in Ref. 3 of the com- 
puted normal force and pitching moment coefficients. 
Such an analysis yields Fourier coefficients that are di- 
rectly related to the normal force and moment stabil- 
ity parameters given by (CW„ - k 2 Cff 4 ), (Cn a + CjvJ, 
(C ma -k 2 C mi ), and (C m& +C m J. The first two terms 
represent pitqh displacement and pitch rate parameters 
while the last two terms represent oscillatory stability 
and damping-in-pitch. Further details of predicting dy- 
namic stability parameters are given in Ref. 3. 

In the current study, the effects of flexibility on 
dynamic stability parameters for the wing-body config- 
uration are examined. To validate the techniques pre- 
sented in this study, a comparison with experimentally 
measured dynamic stability parameters for the rigid 
configuration is given in Fig. 12. In order to match 
with experimental flow conditions, the computations 
are performed at a freestream Mach number of 0.70 
and a pitch amplitude of 1 deg. In Fig. 12, stability 
parameters based on normal force coefficients are given 
for oscillations about various mean angles of attack. 

The effect of flexibility on dynamic stability pa- 
rameters at various Mach numbers are given in Fig. 
13 for mean AOA of 4 deg. Figures 13a- 13d represent 
pitch rate, pitch displacement, damping-in-pitch and 
oscillatory stability parameters, respectively, as defined 
above. It is noted from Fig. 13 that the effects of flex- 
ibility are most pronounced in the pitch displacement 
and oscillatory stability parameters for the entire Mach 
number range. In fact, the effect of in-phase wing flexi- 
bility as defined in this study leads to a decrease in pos- 
itive oscillatory stability {(C m<x - k 2 C m<i ) < 0)* Figure 
13 also shows that more significant changes in stability 
parameters as a function of Mach number occur in the 
transonic regime of M = 0.90 to 1.0. 

Figure 14 illustrates the effect of AoA on dynamic 
stability parameters at a Mach number of 0.95. As in 
Figure 13, more pronounced changes between rigid and 
flexible cases are noted for the pitch displacement and 
oscillatory stability parameters. As in the earlier phase 
and magnitude results, significant non-linear responses 
are noted throughout the angle of attack range. 

Conclusions 

The use of the Euler/Navier- Stokes flow equations 
to compute stability parameters required for design 
are demonstrated. Computations are made at tran- 
sonic/vortical flow conditions which are beyond the lim- 
its of current codes based on the linear aerodynamic 
theory typically used in stability parameter computa- 
tions. Capability of current high fidelity based flow 
codes to account for flexibility is also demonstrated. 

Computations are made for transonic/ vortical 
flows at various Mach numbers and angles of attack for 
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both rigid and flexible cases. It is found that for the 
configuration selected with a highly swept wing, the 
flow characteristics were moderately non-linear with 
respect to Mach number and angle of attack. Even 
at these moderate non-linear flow conditions, the wing 
flexibility influenced the responses and stability param- 
eters considerably. 

Though this work demonstrated the feasibility of 
using the Euler/ Navier-Stokes equations for computing 
the stability parameters, the associated computational 
expense is still an issue. Order of magnitude speed-up is 
required to make this a routine design tool. Such speed- 
up can be achieved by using the advanced computing 
systems such as parallel computers. 
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Fig. 1 The grid topology of the wing-body configuration. 
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Fig. 3 Comparison of displacements at wing tip between 
rigid and flexible cases 
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Fig. 4 Comparison of total lift responses between rigid 
and flexible cases at M = 0.95, a = 0 deg and k = 1. 
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Fig. 2 Comparison between computed and measured 

pressure distributions at M = 0.90 and a ~ 4.09 deg Fig. 5 Comparison of moment responses between rigid 

and flexible cases at M = 0.95, a = 0 deg and k — 1. 







Fig. 6 Effect of Mach number on the magnitude of lift 
for rigid and flexible cases at a = 0 deg and k = 1. 



Fig. 7 Effect of Mach number on the phase angle of lift 
for rigid and flexible cases at a — 0 deg and k =1. 



Fig. 8 Effect of Mach number on the magnitude of lift 
for rigid and flexible cases at a = 4 deg and k = 1 . 



Mach number 


Fig. 9 Effect of Mach number on the phase angle of lift 
for rigid and flexible cases at a = 4 deg and k =1. 



Fig. 10 Effect of angle of attack on the magnitude of 
lift for rigid and flexible cases at M = 0.95 and k — 1. 



Fig. 11 Effect of angle of attack on the phase angle of 
lift for rigid and flexible cases at M = 0.95 and k =1. 
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Pitch displacement parameter Pitch rate parameter 



a, (deg) 

Fig. 12 Comparison of pitch rate and pitch displace- 
ment stability parameters between computations and 
experiment at various angles of attack for M = 0.70 
and k = 0.268. 



Mach number 


Fig. 13-a Effect of Mach number on rigid and flexible 
pitch rate parameter for cases at a = 4 deg and k = 
1 . 0 . 
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Fig. 13-b Effect of Mach number on rigid and flexible 
pitch displacement parameter at a = 4 deg and k = 



Mach number 


Fig. 13-c Effect of Mach number on rigid and flexible 
damping-in-pitch parameter at a = 4 deg and k = 1.0. 



Mach number 

Fig. 13-d Effect of Mach number on rigid and flexible 
oscillatory longitudinal stability parameter at a = 4 
deg and k = 1.0. 
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Fig. 14-a Effect of angle of attack on rigid and flexible 
pitch rate parameter at M = 0.95 and k = 1.0. 



Fig. 14-b Effect of angle of attack on rigid and flexible 
pitch displacement parameter at M = 0.95 and k = 
1 . 0 . 



Fig. 14-c Effect of angle of attack on rigid and flexible 
damping-in pitch parameter at M = 0.95 and k = 1.0. 



Fig. 14 d Effect of angle of attack on rigid and flexible 
oscillatory longitudinal stability parameter at M = 0 95 
and k = 1.0. 
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